THE SCATTERING PROBLEM IN RYCKMAN'S CLASS OF 
JACOBI MATRICES 



L. GOLINSKII, A. KHEIFETS*, P. YUDITSKII** 



Abstract. We give a complete solution of the scattering problem for Jacobi 
matrices from a class which was recently introduced by E. Ryckman. We char- 
acterize the scattering data for this class and illustrate the inverse scattering 
on some simple examples. 



1. Introduction 

In mid 70th Guseinov [8, 9] developed a scattering theory for infinite Jacobi 
matrices 

bi ai 



J = J({a„}, {&„}) = 



a i 




b 2 

« 2 



« 2 



(1.1) 



a n > 0, b n = b n , which can be viewed as a discrete version of the scattering theory 
for one-dimensional Schrddinder operator on the half-line by Marchenko-Faddeev. 
The basic assumption on J is the finiteness of the first moment 



71=1 



n(\a n - 1| + |6 n |) < oo- 



(1.2) 



We say that a Jacobi matrix J = J({a„}, {&„}) belongs to Guseinov's class Q if 
its parameters satisfy (1.2). Later Geronimo [4, 5] (see also [6]) solved the spectral 
problem for Jacobi matrices in more general "weighted" Guseinov's classes by using 
the inverse scattering technique. The main feature of his results is that the decay 
of the Jacobi parameters {a n — 1}, {b n } manifests itself in the decay of the Fourier 
coefficients of the absolutely continuous part (after suitable modifications) of the 
measure. 

In the modern scattering theory of Jacobi operators (see, e.g., [3]) the vari- 
ous analogues of (1.2) (for much more complex backgrounds) appeared, and they 
seemed to be indispensable. 

In 2007 Ryckman [14, 15, 16] came up with a new class of Jacobi matrices, for 
which he obtained a complete spectral description. To state his result we introduce 
some notations and definitions. Let us write 



P = {Me£l S>0 if ||/3| 



< 00. 
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Definition 1.1. A Jacobi matrix J = J({a„}, {b n }) belongs to Ryckman's class 
1Z, or its spectral measure cr(J) G 1Z, if the series ^2 n (a n — 1) and J2 n are 
conditionally summable, and 

oo oo 

■= ~ E bkG ^' Vn — - E - !) G £ ? 

k— n+1 /c— n+1 



Definition 1.2. A function g on the unit circle T is said to be in the Besov class 
B 2 1/2 if the sequence of its Fourier coefficients is in l\ (Z) 

ff(*) = 2sn* n , ^|n|| ff „| 2 <c^. (1.3) 



nGZ 



Also, a function / on the interval [—2, 2] is said to be in -B^ 2 if f(t) f {t + \ ) 
is in -B^ 2 - 

Note that / is a symmetric function, f(t) = f(t), and conversely, each symmetric 
function has the form /. If 

7(*) = E/«*" 

then the symmetry of / implies /_„ = /„. If / is in addition a real function, then 

f—n fn fn- 

Theorem 1.3 (Ryckman). J £ 1Z if and only if the spectral measure <r(J) of J 
has this structure: 

• The absolutely continuous part is supported by [—2,2], and 

a ac ( dx )=f( x ,J)d X ^^lV^^d X , p(x,J)= i2 _ P x P { zl xp2 (1-4) 

l li 

with 71,72 equal or 1, and logpo eS 2 ' . 

• The singular part is 

N 

a s (dx)=J2^k(J)S(X k ), N = N(J)<oo, <i k (J) > 0, A fe e R\[-2, 2]. (1.5) 
fe=i 

Note that G ClZ, and the inclusion is proper. Indeed, 

oo oo / oo \ ^ oo oo oo 

5>k™i 2 =X> E M =EEE«nn 

n—1 n—l \k— n+1 / n— 1 fc= n /— n 

oo oo oo oo oo oo 

^EEE z i^ii 6 ^EEE z i^ii 6 'i 

n—l fc— n l—n n—l k—n l — l 

oo k oo oo oo / oo 

=EEE z NM = E fc E«<i= E fc N 

fc=ln=i;=l fe=l (=1 \fe=l 
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Similarly, 



5>M 2 < f 5^As|o fc — 1| 

71=1 \fe=l 
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On the other hand, J{{a n },{b n }) with 

a n = l + (-l) n n- 1 - e , b n = {-l) n n- 1 ~\ 0<e<l 

belongs to 1Z, but (1.2) is false. 

Given a(J) £ 1Z, define two outer functions which are key ingredients in both 
spectral and scattering theory of the class 1Z. First, we put 



D (z) =D (z,J) :=cxp|i J log p {t,J)m(dt)^ 
u (z) + iv {z) ' 



(1.6) 

] u (z) -+- iv \z) i 
= CXP { 2 } 

with 

uq{z)= I P z (t)logp (t, J)m(dt), v (z) = / Q z (t) log p (t, J)m(dt) = u (z) 
Jt jt 

a harmonic conjugate to uq, 

u (t) = logpo(t, J) = ^2r k t k , v (t) = ^^2(sgn k)r k t k , 

fcez ' fcez 

both Uq and Vq are real, mq is symmetric, and vq antisymmetric: Vo(t) = —v (t). So 



D (t) = D (t), |Do(*)r = Pb(*> J ) (1-7) 



almost everywhere on T. It is known (cf., e.g., [17, Proposition 6.1.5]) that log/So € 
^ 2 implies (po) ±1 
Secondly, we put 



B 2 1/2 implies (p ) ±1 € £ P (T) for p < oo, so D* 1 e i? p for such p. 



£>(*) - £>(*, J) := exp |1 jf logp(i, J)m(dt) J 



Do(z,J) (L8) 



(1 - z) 7l (l + z)T2 ' 



Both Dq and D are related to the absolutely continuous part of the spectral 
measure. The discrete part is completely determined by the set of eigenvalues 
{Xk}, or equivalently, by the set 

Z{J):= \z k (J):X k =z k (J) + ^—, fc = l,2,...,Jvl ) (1.9) 

z k (J) £ (—1, 1)\{0}, and by the set of masses {cr k {J)} k=1 in (1.5). 

Definition 1.4. Given J £ 1Z, under the scattering data for J we mean the follow- 
ing collection {71 ( J), j 2 ( J); Z(J); pi(J), . . . , p N (J); s(t,J)} 

(1) A pair ( 7l (J), 72 (J)) from (1.4); 

(2) The set Z(J) from (1.9), or equivalently, a finite Blaschke product 

**^-n £(V-~»(% - ^"^^^Blbj- teT; <L10) 

(3) N = N(J) positive numbers 



p k {J) := <r k {J) 



B'{z k ,J) 
D(z k ,J) 
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\l-z k {J)- 2 \ >0, k = l,2,...,N; (1.11) 
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(4) The scattering function 

0o is the Jost function for J (see Section 2). 

Compared to [4, 5] we move backward, from spectral to scattering. The goal of 
the present note is to obtain a complete characterization of the scattering data in 
Ryckman's class, and so demonstrate that the scattering theory goes far beyond 
Guseinov's class (1.2). We analyze the scattering data, prove the uniqueness the- 
orem in Section 2, and solve the inverse scattering problem in Section 3. A few 
simple examples are given in Section 4. 

2. Scattering data 

The basic three-term recurrence relation for a Jacobi matrix J({a n }, {b n }) 

a„_i2/„_i + b n y n + a n y n+1 = (z + z~ x )y n , n = l,2, a = 1 

has two "distinguished" solutions. The first one, known as the sine-type solution, 
is 

Vn = S n {z) = p n -i I Z + - ) , S =0, Si = 1, 



Z, 

Pk are orthonormal polynomials with respect to the spectral measure cr(J). A 
fundamental result by Szego concerns an asymptotic behavior of orthonormal poly- 
nomials with respect to "nice" measures with supper C [— 2,2]. It was extended 
substantially in [7, 11] and [12], where a finite (respectively, infinite) number of 
mass points outside [—2,2] is allowed. In our notation the Szego asymptotics states 
that for J e K 

Q(z) := lim z n p n (z+-) = - B< £$ n (2.1) 

n^oo \ z) (1 — Z Z )D(Z, J) 

uniformly on the compact subsets of the unit disk D. 

The second solution is the Jost solution y n — <j> n , defined by a specific asymptotic 
behavior at infinity 

lim z^ n (j) n (z 1 J) = 1 

n— >oo 

uniformly on the compact subsets of the unit disk D. The Jost solution exists under 
certain additional assumptions (cf. [18, formulae (13.9.2)-(13.9.4)]), which are met 
for J e 1Z. <fio is called the Jost function. For an exhaustive treatment of the Szego 
and the Jost asymptotics see [1], [18, Section 13.9]. 

The relation between the Szego asymptotics and the Jost function is given by 
(see, e.g., [2], [18, Theorem 13.9.2])) 

M z) = (l-z*)Q{z) = ^l. (2.2) 

Hence, for J e 1Z (f>o € H p with p < oo, so the boundary values of 0o exist a.e., 
and the scattering function s (1.12) is well defined. Moreover, by (1.8) 



_ (i-t)^>(i+*r( j > p (i,j) 2 

l ' > (l-t)Ti(J)(l + *)'»(■') D (t,J) l ' > 

D (t, J) 



(2.3) 
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Clearly, s(t) = s(t) = s" 1 ^). 

1 /2 

Theorem 2.1. The scattering function of a Jacobi matrix J <G 1Z belongs to B 2 ' 
and admits representation 

s{t,J) = (-l)^ J h M e- lvW , (2.4) 
where M = 2N + 71 (J) + 72 (J) G Z +) u satisfies 

v(t) = : ^At) = -v(F), uefif. (2.5) 
Proof. Wc apply (2.3). By (1.6) 



goM) = = IA)(M)| 2 = e -i«o(t) 

D (t,J) D {t, J) D 2 Q (t,J) 

Since uo = logpo € B^ 2 ' an d the Hilbert transform is bounded (isometric) in B^ 2 , 
then vq e B^ 2 and antisymmetric. 
Next, 

B 2 (t,J)=t 2N (f[l-^-) =t ™e- i ^ t \ 

Vfe=i 1 Zfc v 

1 /2 

and as Zfe € D then v\ e £? 2 ' and antisymmetric. It remains to put w = vq + v±. 

1/2 

It is clear from (2.4) that s € B 2 , as claimed. □ 

Remark 2.2. Representation (2.4) for the scattering function of J <G 7\L is a direct 

^ 1/2 

consequence of (2.3) and logpo & B 2 . A slight refinement of Peller's theorem [13, 
Corollary 7.8.2] states that an arbitrary function 

heBl /2 , h(t) = h{Fj = h- 1 (t) 

a.e. admits representation 

h(t) = (-iyt ] e- tw ^\ 

where 7 = or 1, j £ Z an integer number, w satisfies (2.5), and such representation 
is unique. A pair (7, j) can be viewed as an index of h. 

Let us turn to the numbers /J, k (J) (1-11). In his version of the scattering theory 
for Jacobi matrices (1.2) Guseinov suggested the normalizing constants 



m k 



(J) :=]T|0„(z/f(J),J)| 2 , fc = l,2,...,JV (2.6) 



as a part of the scattering data. We show that these values agree. 

Proposition 2.3. LetJeK. Then fi k (J) = m k (J), k = 1, 2, . . . , N. 

Proof. It is known from the general theory of Jacobi matrices and orthogonal poly- 
nomials, that the vectors 

life = {s n (z k {J))}„>i = {p„(Afe)}„> € £ 2 , 

so life are eigenvectors of J with the corresponding eigenvalues Afe. Furthermore, 

1 00 00 

-JJ-^Yl \Sn(MJ))\ 2 = ]T \ Pn (\ k )\ 2 . (2.7) 
ak[ - J > n=l n=0 



6 



On the other hand, 4>a(z k (J)) = 0, and so $ k = {<fi n (zk(J), J)}n>i are also eigen- 
vectors of J for the same eigenvalues. Hence = Cfcllfc, and we find the constants 
Cfe from the initial data Si = 1, so that c k = <fii(z k (J), J)- By (2.6) and (2.7) 

mk(J) = |&(* fc (J),.7)|W)- 

It remains to express </>i in terms of the spectral data. Once the Jost asymptotics 
exists for J € 1Z, the Jost solution </>„ is proportional to the Weyl solution 

w n {z) ■= ((z + z~ x - J) _1 ei,e„) , n=l,2,..., io = 1, 

that is, </>„ = <pow n . In particular, 

^i(z, J) = <^o(2, J)wi(z) = o (z, J)M(z, J), 

where M is the Weyl function for J 

m ( >, j) = + fl-' ei , ei ) . / r . _fiW_ + sr (l) , 

M is analytic at z k {J). So 

Since lim z _ >Zfc (z — z k (J))M(z, J) = a k (J)(l — z / 7 2 (J)) _1 , we finally have 

<pl{Z k {J), J) = , 

1 - z fe (J) L>{z k {J),J) 
as needed. □ 

To complete the analysis of scattering data we prove the uniqueness theorem. 
Theorem 2.4. Let Ji £ 1Z, I = 1, 2, Ziaue £/ie same scattering data. Then J\ = J 2 . 
Proof. We want to make sure that cr(Ji) = o-(J 2 ). It is clear from (2.3) that 



Dpjt, Ji) = D (t,J 2 ) 
A, (Mi) D (t,J 2 )' 
But in Ryckman's class Dq 1 e H 2 , so the latter means D (J 2 ) = cDq(Ji), D(J 2 ) = 
cD(Ji) for some c > 0, and hence a ac (J2) = c 2 a ac (Ji). Next, fi k (Ji) — HkiJi) 
implies by (1.11) cr k (J 2 ) = c 2 cr k (Ji), and the normalizing condition 
,2 n 2 at 

/ f{x,J 1 )dx + '^2a k (J 1 )= f(x, J 2 ) dx + y~] <7 fc (J 2 ) = 1 
•^~ 2 fe=i fe=i 

gives c = 1, as needed. □ 

3. Inverse scattering 

Consider the following collection of data {71, 72; Z; fi^', s}' 

(1) a pair of numbers (71,72) from {0, 1} x {0, 1}; 

(2) an arbitrary set of N distinct points Z = {z k } k=1 in (—1, 1)\{0}; 

(3) an arbitrary set of N positive numbers /i k ; 

1/2 

(4) a function s <G B 2 , \s\ = 1 a.e. on T, with the index (71, 2N + 71 + 72), 
i.e., 

s(t) = (-1)^ t 2N+ ~> l+ ~* 2 e~ l " (t) , 
uj satisfies u){t) = u){t) = —u(t), lj € B^ 2 . 
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Theorem 3.1. There exists a unique Jacobi matrix J € 1Z, for which the above 
collection is the scattering data. 

Proof. As in the proof of Theorem 2.1 we can write 

( ) (1-4)71(1 + a [t ^ )e 
vq is subject to (2.5). The Fourier series for vq is 

Vo(t) = *^2v (n)t n , vo(-n) = v (n) = -v (n), 



so ub(0) = 0. Take uq such that vo is its harmonic conjugate. Then Uo(t) = uo(t) = 
u (t) and 

uo(t) = u (n)t n , uo(-n) = u (n) = v (n). 

Note that uq is defined up to an additive real constant Mo(0), which will be chosen 
later on from the normalization condition. 

Define a function p on [—2, 2] by p = e u ° , and put 

: = • m == s*) 

both up to a factor C = e M °(°). Next, write 



r> ( \ (I [t + Z \ ( U Q (z)+iv Q (z)\ 

A)0)=expj- / j— u a {t)m{dt)\ = expj J, 



(1-Z)T1(1 + Z)- 

and put 



o-fe := Hk 



D{z k ) 



2 



1-Zfe 2 | ->0, k = l,2,...,N, 



-2 



B'(z k ) 

the latter values are defined up to a factor C above, which is now taken from 

N 



/ f{x)dx + y^a k = 1. 
J ~ 2 fc=l 



1 /2 

Since i>o € -B 2 j then so is uo, and by Ryckman's theorem the measure cr = 
{/, {<7fc}} is the spectral measure of some Jacobi matrix J E 1Z. By construction, 
{71, 72; Z; hi, . . . , hn', s} is the scattering data for J, and J is unique by Theorem 
2.4. The proof is complete. □ 

4. Examples 

As we see, it is comparatively easy to restore the spectral measure from the 
scattering data. Assume first that N = 0, that is, supper C [—2,2]. To find the 
Jacobi parameters it seems reasonable now to carry the measure over from [—2,2] 
to the unit circle (the inverse Szego transform) a = Sz^°\h) 1 

o-(dx) = f{x) dx, jj,{dt) = w(t)m(dt), w(t) = , (4.1) 



1 Due to the form of (1.4) it is convenient to use the modified Szego transform 
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(i is a probability measure on T (by c we denote different positive normalizing 
constants), to compute the Verblunsky coefficients a n — a n {p) of fi, and then go 
back to Jacobi parameters with the help of the Geronimus relations 

bn+l =Q!2n(l - Ck2n+l) - a2n+2^ + a 2n+ l), 

(4.2) 

a 2 n+ i={l-a 2n+ z{l-al n+2 )(l + a 2n+ i), n = 0, 1, 
see [18, Section 13.2]. 

Example 4.1. Let 71 = 72 = 0, a e [0, 1), and s(i) = (1 — at)(l - at) -1 . Since 
B = lwe have as in the proof of Theorem 3.1 



r 2 



1 — az ' |1 — at\ 2 1 — ax + a 2 

so by (4.1) 

/~"m ( rii~ 1 



|l-at| 2 ' 

Hence is the Bernstein-Szego measure, for which the Verblunsky coefficients are 

ao = a, cti = a 2 = • • • = 0. 

Finally, by (4.2) 

61 = a, b 2 = 63 = • ■ • = 0, ai = a 2 = . . . = 1. 

Example 4.2. For 71, 72, a as above put s(t) = (1 — ai)(l — at) -1 . We now have 

£>(*) = c(l - az) , p(t) = c|l - at| 2 , f(x) = c(l - ax + a 2 )V4 - x 2 , 

and //(dt) = cj 1 — at| 2 m(dt). The Verblunsky coefficients are (see [17, Example 
1.6.4]) 

a 1 — a 

«n = ~ a _ n _ 2 _ an+2 , n = 0,l,.... 

Finally, 

b , , - n 2n+1 ^ - o 2 - 1 „ 2 "+ 2 Z q2 ) 2 

° (1 _ a 2 n +2) (1 _ a 2n+4) ' a n+l-± « (1 _ a 2„+4)2 ' 



It is worth pointing at the difference between the above examples. In the first one 
D^ 1 is a polynomial (of degree 1), so there are finitely many nonzero Verblunsky 
coefficients a n , and J is of finite support, i.e., a n = 1 and b n = for all large 
enough n. In the second case D~ Y is a rational function with the pole 1/a, so the 
Jacobi parameters tend to their limits exponentially fast. 

Example 4.3. Let 71 = 72 = N = 0, a, b e [0, 1). Put 

(l-at)(l-bt) 
W (1 - at)(l - bt) ' 
As above, the spectral measure is 



. . cV 4 — x 2 , , . cm(dt) 

<7(dx) = — — — — dx, n(dt) 



{l-ax + a 2 ){l-bx + b 2 ) ' ^ y ' |l-at| 2 |l-6t| 2 ' 
The latter is again the Bernstein-Szego measure with the Verblunsky coefficients 
a + b 

a = — — - , ai = -ab, a 2 = a 3 = . . . = 0. 
1 + ao 
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Finally, 

bi — a + b, 62 = ^3 = • • • = 0; a 2 = 1 — ab, a 2 = as = . . . = 1. 

Remark 4.4. For 71 = 72 = the Szego transform Sz^ -* is a right one in the sense 
of [14]. We can compute the same examples with 71,72 = 0, 1 by using all four 
Szego transforms Sz^°\ Sz^°\ Sz^, correspondingly. 

To complete the examples section consider the case N = 1, that is, the spectral 
measure has a mass point outside [—2, 2]. Now the Szego transforms do not work, 
so we proceed in two steps. First, we compute the Jacobi parameters a n ,b n for 
<j\ = ca ac (properly normalized), and then we add a mass point and recompute the 
Jacobi parameters (cf., e.g., [10, Lemma 7.15]). 

Example 4.5. Let 71 = 72 = 0, N = 1, z\ € (0,1) and \i\ > are given, and 
s(t) = t 2 . We have 

j. \ 2 / 1 T \ 2 

2/ \ / — 4-2 1 — 1 

so 



B (t) = — - ) =t ■ 



'D(t)' w (1-zxz) 2 ' 

The spectral measure is of the form 



/ , x Cqa/4 - x 2 , , 1 

g(d!c) = 27r(l-z 1 x + z 1 ^ fr + ^fr), = + - 

with 



(Tl = 111 



D( Zl ) 



|i-*rT = cg£ 



B'(zO 

c is determined from cr(R) = 1, or 

f2 , , , , Lll ... \/4 — X 2 



fo(x)dx + -4, / (a;) 
-2 z i 



2tt(1 - zirr + z 2 ) 2 ' 

Take oo(<ix) = effodx, the measure on [—2, 2], ero(R) = 1- It is not hard to find 
the values of normalizing constants co,ci: 

'2 V4~ 1? 

— — ax. = — 

.2 



_J_ Mi a / 2 dx 1 



If Hi(dt) = wi(t)m(dt) and (T = Sz^(//i), then w\ = c|l — zii| 4 , and as in 
Example 4.3 (with a = b = z\) 

&i(<ro) = 2*i, & 2 (<7 ) = ... = 0; a 2 (<7 ) - 1 - z 2 , a 2 (<r ) = . . . = 1. (4.3) 

Note that <7o is the Bernstein-Szego measure on [—2,2] (see [19, Section II. 2. 6]), 
and the corresponding orthonormal polynomials are known explicitly 

z n+1 (l- - z^-^l- Zl z) 2 
p n (a a ,z + z ) = z , n=l,2,..., p = 1. 



In particular, pk{<Jn, Ai) = ^/l — z 2 z 1 , fc = 1,2, The Christoffel kernels are 

n 

tf„+i(<7 ,Ai) =5Zpfe(<7 ,Ai) = zf 2 ™, n = 0,l,..., ^o = 0. 



fc=0 
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so 



To apply Nevai's formulae write the perturbed measure a in a canonical form 

a(d X ) = ^ dX ]l £S ^\ e = ff- 1 = ^(1-,?), (4.4) 
1 -\- s Cq z^ 

2/ \ 2/ \ (! +e-f^n-i(o-o, Al))(l +£-K n+ i(o- , Ai)) 



(1 + AO) 2 
M ") =6„(ct ) - an-i(ffo)K-i + an(o-o)K, n = 1,2,. 



where 



)£(!— ^i) \ 9 

Zl(z ^-2 +£) , « ^ ^; _ 
n-l ' 
zi(l+e) ' " ±- 



Eventually, we have 



,. 2 l+Zi^l + Z! 2 ) 

(i + Ml zr 4 (i-^)) 5 



1 _1_ £ (1 Z l) 2 ,2n ^ r, o 
Q r» =1+ / 2n , g ^N2 Z l ' ™ = 2 , 3, • • • 



and 



i, /„\ £ (1 z l) 2 „2n-5 \, q 



(e + zl n - 2 ){e + z?«- 4 ) 1 

£ from (4.4). 

Remark 4.6. The Jost function is now 

(f>o{z) = c(l - ziz) ( 1 - — 
V z i 

but J is of infinite support. The Jacobi parameters tend to their limits exponentially 
fast, cf. [2, Remark 1.10]. 
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